N  PS  ARCHIVE 
1969 
DESPARD,  R. 


LAMINAR  BOUNDARY   LAYER  SEPARATION    IN 
OSCILLATING  FLOW 

by 

Ronald  Arthur  Despard 


DUDLEY  KNOX  LIBRARY 

NAVAL  POSTGRADUATE  SCHOOL 

MONTEREY,  CA  93943-5101 


United    States 
Naval  Postgraduate   School 


THESIS 


LAMINAR  BOUNDARY  LAYER  SEPARATION 
IN  OSCILLATING  FLOW 


by 


Ronald  Arthur  Despard 


June  1969 


TIiua  documznX.  ka6  bzzn  appiovzd  ^ok  pubtic  n.z- 
Izoaz  and  6 ale.;  -LU  dl&VUbiUion  -ca  unlonitzd. 


i"s1.'al/avalPostSriduate9Sc»oo1 


UontexejL. 


California 


Laminar  Boundary  Layer  Separation 
in  Oscillating  Flow 


by 


Ronald  Arthur  Despard 
Lieutenant  Commander,  United  States  Navy 
B.A.E.,  Rensselaer  Polytechnic  Institute,  1959 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


DOCTOR  OF  PHILOSOPHY 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
June  19  6  9 


ABSTRACT 


Laminar  boundary  layer  separation  phenomena  in  oscil- 
lating flow  was  experimentally  investigated.   Multiple  hot- 
wire anemometers  were  used  to  obtain  instantaneous  boundary 
layer  velocity  profiles  on  a  model  in  an  oscillating 
freestream.   Certain  instantaneous  profile  behavior  was 
found  to  be  uniquely  related  to  wake  formation,  while  non- 
wake-inducing  transient  flow  reversals  were  found  to  occur 
throughout  the  adverse  pressure  gradient  regime.   Based  on 
these  observations,  a  practical  definition  of  boundary 
layer  separation  in  oscillating  flows  was  formulated. 
Separation  data  then  obtained  revealed  that  the  imposition 
of  oscillations  caused  the  separation  point  to  move  upstream 
of  its  steady  flow  location.   From  the  data  acquired,  a 
parametric  correlation  of  separation  point  location  was 
developed. 


Llbrary  DUDLEY  KNOX  LIBRARY 

US.   Naval  Postgraduate   q0u  naVAL  POSTGRADUATE  SCHOOL 
Monterey,    Calif o^  ^J^TEREY,  CA  93943-5101 

TABLE  OF  CONTENTS 

I.  UNSTEADY  LAMINAR  BOUNDARY  LAYER  FLOWS  17 

A.  INTRODUCTION 17 

B.  PREVIOUS  ANALYTICAL  STUDIES  17 

C.  PREVIOUS  EXPERIMENTAL  STUDIES  20 

D.  TRANSITION  AND  SEPARATION 22 

II.  LAMINAR  BOUNDARY  LAYER  SEPARATION 

IN  OSCILLATING  FLOWS  26 

A.  INTRODUCTION 26 

B.  CONCERNING  THE  DEFINITION  OF  SEPARATION  26 

C.  OBJECTIVES 28 

III.  EXPERIMENTAL  EQUIPMENT  29 

A.  WIND  TUNNEL 29 

1.  General  Description  29 

2.  Rotating  Shutter  Valve — -  31 

3.  Test  Section 32 

B.  MODEL 35 

C.  INSTRUMENTATION 40 

1.  Freestream  Sensors "■"""•■ ■■""- — .........  -43 

2.  Hot-Wire  Anemometers  40 

D.  MULTIPLE  HOT-WIRE  PROBE  45 

1.  General  Characteristics 45 

2.  Probe  Modules  48 

3.  Vertical  Traverse  Probe  50 

4.  Probe  Carriage 50 


IV.  EXPERIMENTAL  PROCEDURE  59 

A.  CHARACTER  OF  THE  FLOW  ENVIRONMENT 59 

B.  INVESTIGATION  PROCEDURE  60 

C.  FLOW  VARIABLES 61 

1.  Mean  Velocity 61 

2.  Static  Pressure  Distribution  61 

3.  Frequency  of  Oscillation 63 

4.  Amplitude  of  Oscillation 64 

D.  DIMENSIONLESS  FLOW  PARAMETERS  64 

V.  RESULTS 66 

A.  FLOW  PHENOMENA 66 

1.  Flow  Region  Classification 66 

2.  Flow  Reversal 68 

3.  Periodic  Effects  71 

B.  INSTANTANEOUS  VELOCITY  PROFILES  72 

C.  DEFINITION  OF  SEPARATION 85 

D.  CORRELATION  OF  SEPARATION  POINT 

ALTERATION  106 

VI.  DISCUSSION 118 

A.  PERIODIC  WAVEFORM  DISTORTION   118 

1.  Amplification 118 

2.  Phase  Shifts 120 

B.  TRANSIENT  FLOW  REVERSAL  122 

VII.  CONCLUSIONS  124 

APPENDIX  A  -  Representative  Profile  Data  125 

APPENDIX  B  -  Mathematical  Outline  of 

Parametric  Correlation 153 


LIST  OF  REFERENCES 157 

INITIAL  DISTRIBUTION  LIST  160 

FORM  DD  1473 161 


LIST  OF  TABLES 


TABLE  PAGE 

I.   Operating  Conditions  for  Instantaneous 

Velocity  Profile  Runs  ■ 73 

II.   Steady  Flow  Separation  Data 107 

III.   Operating  Conditions  and  Flow  Parameters 

at  Separation 108 


LIST  OF  ILLUSTRATIONS 

FIGURE  PAGE 

1.  Plan  View  of  the  Wind  Tunnel 30 

2.  Photograph  of  the  Rotating  Shutter  Valve  32 

3.  Photograph  of  the  Test  Section 34 

4.  Typical  Wind  Tunnel  Test  Section 

Velocity  Profile  36 

5.  Overall  Photographic  View  of  the 

Wind  Tunnel : 37 

6.  Photograph  of  the  Model 38 

7.  Photograph  of  the  Ten-Channel 

Hot-Wire  Anemometer  41 

8.  Typical  Hot-Wire  Anemometer  Calibration 

Curve 43 

9.  Typical  Hot-Wire  Anemometer  Calibration 

Velocity  Profile  in  Blasius  Flow  44 

10.  Schematic  Diagram  of  Instrumentation  46 

11.  Photograph  of  the  Control  Console 47 

12.  Photograph  of  a  Probe  Module 49 

13.  Photograph  of  the  Vertical  Traverse 

Probe  Disassembled 51 

14.  Photograph  of  the  Vertical  Traverse 

Micrometer  Control  52 

15.  Photograph  of  the  Probe  Carriage  with 

Modules  and  Vertical  Traverse 

Probe  Mounted 53 

16.  Effect  of  Presence  of  Probe  Carriage 

on  Velocity  Profile  55 

17.  Photograph  of  the  Chordwise  Traverse 

Mechanism  Mounted  in  the  Test 

Section 57 


18.  Photograph  of  the  Probe  Assembly  and 

Model  Mounted  in  the  Test  Section 58 

19.  Mean  Pressure  Coefficient  Distribution  62 

20.  Typical  Velocity  Waveforms  in  the  Three 

Adverse  Pressure  Gradient  Flow 

Regions 67 

21.  Oscillating  Velocity  Component  Waveforms  70 

22.  Typical  Sequence  of  Instantaneous  Profiles 

during  Half-Cycle  of  Oscillation, 

Pre-Wake  Region  (Weak  Transient 

Flow  Reversal)  77 

23.  Oscillating  Velocity  Component  Waveforms, 

Pre-Wake  Region  78 

24.  Typical  Sequence  of  Instantaneous  Profiles 

during  Half-Cycle  of  Oscillation, 
Pre-Wake  Region  (Strong  Transient 
Flow  Reversal)  80 

25.  Oscillating  Velocity  Component  Waveforms, 

Pre-Wake  Region  81 

26.  Examples  of  Streamwise  Development  of 

Reverse  Flow  in  the  Cycle  of 

Oscillation 83 

27.  Typical  Variation  of  Instantaneous  Profile 

Envelope  with  Increasing  Reynolds  Number 84 

28.  Typical  Cycle  of  Instantaneous  Velocity 

Profiles  at  Separation  (Weak 

Flow  Reversal)  88 

29.  Oscillating  Velocity  Component  Waveforms 

at  Separation 89 

30.  Typical  Cycle  of  Instantaneous  Velocity 

Profiles  at  Separation  (Strong 

Flow  Reversal)  90 

31.  Oscillating  Velocity  Component  Waveforms 

at  Separation 91 

32-58.  Envelopes  of  Instantaneous  Profiles  at 
Separation  (for  various  operating 
conditions)  92 


10 


59-62.  Separation  Point  Alteration  as  a  Function 
of  Pressure  Parameter  and  Reynolds 
Number  (for  values  of  frequency 
parameter)  114 

63.  Variation  of  Maximum  Amplification 

with  NA 119 

64.  Variation  of  Maximum  Phase  Lead 

with  N,,  121 

A 


11 


Symbol 

c 

S 

E 
u 

NA 

Nf 

N„ 

LIST  OF  SYMBOLS 

Definition  Units 

Chord  length  inches 

Mean  static  pressure  coefficient  

Electrical  signal  corresponding  to  U  volts 

Dimensionless  amplitude  =  U/U  


Dimensionless  frequency  parameter  =  — — =- 

U  Z 
.  °° 

Dimensionless  pressure  gradient- 
history  parameter 


x  _ 

J  x  -s—   dx 

x_  dx 
o 


/    x 


dP 


dx 


dP 


dx 


x 

o 


■Jx    q<x> 

dx      s 


N      Length  Reynolds  Number 
Ke 


u  2 


U  X 

oo 

V 

P      Mean  static  pressure  lb. /in 

P      Remote- free stream  mean  static  pressure  lb. /in 


Remote-f reestream  mean  dynamic  pressure     lb. /in 

=  P 


2tt 

T      Period  of  oscillation  =  —  sec. 

(0 

t      Time  measured  from  maxima  in  sec. 
freestream  velocity 

U      Instantaneous  local  freestream  velocity  ft. /sec. 

U^     Instantaneous  remote  freestream  velocity  ft. /sec. 

U      Local  mean  freestream  velocity  ft. /sec. 

U      Oscillating  freestream  velocity  component  ft. /sec. 

U      Remote  mean  freestream  velocity  ft. /sec. 


13 


u      Instantaneous  boundary  layer  velocity      ft. /sec. 

u      Mean  boundary  layer  velocity  ft. /sec. 

u      Oscillating  boundary  layer  ft. /sec. 

velocity  component 

x      Distance  from  leading  edge  in. 

x      Point  at  which  the  surface  pressure  in. 

gradient  becomes  positive  for  the 
last  time  prior  to  separation 

x      Location  of  oscillating  flow  in. 

s 

separation  point 

x      Location  of  steady  flow  separation  in. 

point 

y      Distance  from  the  wall  in. 

A      Dimensionless  alteration  of  separation     

s 

point  location 

x   -  x 

_   ss s 

xss 
6       99%  boundary  layer  thickness                 in. 
ri      Blasius  parameter  


<J>      Phase  angle  degrees 

2 
v      Kinematic  viscosity  ft. /sec. 

3 
p      Density  lb. /ft. 

a)      Circular  frequency  1/sec. 


14 


ACKNOWLEDGEMENTS 

The  author  wishes  to  express  his  sincere  gratitude  to 
Dr.  James  A.  Miller,  friend  and  adviser,  for  his  interest, 
help  and  guidance  in  the  preparation  of  this  work. 

Grateful  acknowledgement  is  also  due  the  entire 
technical  engineering  staff  of  the  Department  of  Aero- 
nautics, and  in  particular  Messrs.  Theodore  B.  Dunton, 
Robert  A.  Besel,  Bertis  Funk,  and  Norman  Leckenby,  for 
their  constant  cooperation  and  assistance  in  this  investi- 
gation. 

Finally,  the  author  wishes  to  thank  his  wife,  Gretchen, 
whose  patience,  help,  and  encouragement  were  major  factors 
in  the  completion  of  this  work. 


15 


I.   UNSTEADY  LAMINAR  BOUNDARY  LAYER  FLOWS 

A.  INTRODUCTION 

Knowledge  of  the  character  and  behavior  of  the  boundary 
layer  in  unsteady  flow  is  of  prime  importance  to  the  full 
understanding  of  many  practical  problems  in  fluid  mechanics. 
Historically,  rotating  stall  in  turbo-machinery,  aero- 
dynamic flutter,  and  the  flow  around  accelerating  bodies  have 
been  the  most  prominent  of  the  problems  in  this  class. 
Recent  increases  in  use  and  rapid  development  of  rotary  wing 
aircraft,  with  their  inherently  complex  aerodynamics,  have 
re-emphasized  the  need  for  such  knowledge. 

B.  PREVIOUS  ANALYTICAL  STUDIES 

Most  of  the  significant  contributions  to  the  resolution 
of  these  problems  are  analytical  investigations  into  the 
response  of  boundary  layers  to  large  scale  freestream 
unsteadiness  resulting  from  accelerating  or  periodically 
fluctuating  flows.   The  earliest  analysis  of  such  flows 
is  the  classic  work  of  Stokes  [1]  who  treated  an  infinite 
lamina  oscillating  in  its  own  plane  in  a  fluid  at  rest. 
Rayleigh  [2]  considered  the  second  order  effects  in  this 
same  problem.   Schlichting  [3]  expanded  Rayleigh' s  analysis 
by  employing  boundary  layer  assumptions.   These  early 
analyses  treated  what  is,  in  essence,  the  simplest  special 
case  of  the  general  problem,  unsteady  viscous  flow  in  the 
absence  of  a  mean  flow  or  pressure  gradients. 
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The  first  significant  analysis  of  unsteady  boundary 
layer  flows  that  included  the  presence  of  a  mean  flow  was 
that  of  Lighthill  [4]  who  treated  the  case  of  small,  low 
frequency  sinusoidal  oscillations  superimposed  on  the  mean 
flows  about  a  cylindrical  body,  and  a  semi-infinite  flat 
plate.   His  analysis  was  based  on  a  "small  perturbation" 
treatment  of  the  boundary  layer  equations  for  sinusoidal 
flow,  retaining  only  first  order  perturbation  terms.   For 
the  lower  ranges  of  frequency,  a  solution  was  obtained  using 
a  Karman-Pohlhausen  method.   For  higher  frequencies,  the 
relatively  small  convective  terms  were  ignored,  resulting 
in  equations  similar  to  those  of  Stokes'  "zero  mean  flow" 
solution.   A  number  of  fundamentally  important  results  were 
determined:   at  low  frequencies,  the  boundary  layer  flow 
is  essentially  quasi-steady,  having  the  same  character- 
istics for  any  instantaneous  magnitude  of  freestream 
velocity  as  it  would  have  in  a  steady  flow  of  the  same 
velocity,  and  is  independent  of  the  dynamic  history  of  the 
freestream  flow.   As  the  frequency  of  oscillation  is 
increased,  a  limiting  frequency  is  reached,  beyond  which 
the  boundary  layer  responds  to  the  oscillations  as  it 
would  in  the  absence  of  a  mean  flow.   This  limiting  solution 
is  commonly  known  as  the  "shear  wave"  solution.   Lighthill 
found  that  the  maximum  value  of  skin  friction  at  any  point 
leads  the  maximum  in  freestream  velocity,  presumably  due  to 
a  faster  response  of  the  lower  velocity  laminae  in  the 
boundary  layer  to  changes  in  the  time-dependent  pressure 
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gradient  than  the  response  of  the  freestream  to  these 
changes.   Lighthill's  work  has  been  treated  as  a  point  of 
departure  for  many  subsequent  investigations.   Lin  [5]  has 
analyzed  the  high  frequency  case  for  large  amplitude 
oscillations  and  confirmed  the  existence  of  the  shear 
wave  solution.   He  also  concluded  that  the  governing 
equations  become  essentially  linear  at  high  frequency, 
allowing  separation  of  the  effects  of  the  oscillatory  and 
steady  components  of  the  flow.   Carrier  and  Di  Prima  [6], 
using  an  analysis  based  on  a  system  of  equations  derived 
by  a  modified  Oseen  linearization  of  the  equations  of 
motion,  confirmed  both  the  shear  wave  solution  and  the 
phase  advance  of  wall  shear  found  by  Ligh thill.   Nickerson 
[7]  has  expanded  the  scope  of  Lighthill's  flat  plate 
analysis  by  initially  retaining  perturbation  terms  beyond 
first  order.   Quite  recently,  Barriol  and  Lucius  [8]  have 
used  numerical  methods  to  obtain  asymptotic  solutions  to 
the  boundary  layer  equations  for  oscillating  flow,  remote 
from  the  leading  edge  of  a  semi-infinite  flat  plate,  in 
the  absence  of  a  pressure  gradient.   Their  solution  is  in 
close  agreement  with  the  results  of  Lighthill  and  Nickerson, 
but  has  a  much  broader  range  of  applicability. 

The  response  of  the  laminar  boundary  layer  to  small 
freestream  oscillations  in  the  vicinity  of  a  stagnation 
point  has  been  investigated  by  both  Rott  [9]  and  Glauert 
[10].   The  exact  solution  for  a  stagnation  flow  is 
demonstrated  by  Glauert  to  be  obtainable  from  a  single 
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ordinary  differential  equation.   Glauert's  solutions  for 
high  and  low  frequency  limiting  cases,  as  well  as  skin 
friction  phase  angle  lead,  agree  with  Lighthill's 
calculations . 

Compressibility  effects  in  the  laminar  boundary  layer 
for  nearly  quasi-steady  flow  were  considered  by  Ostrach 
[11] ,  who  treated  first  order  deviations  from  quasi-steady 
characteristics.   His  work  again  confirmed  the  skin 
friction  phase  angle  lead  found  by  others.   Moore  [12] 
improved  the  analysis  by  including  higher  order  effects 
for  both  isothermal  and  adiabatic  flat  plates.   He  obtained 
time-average  velocity  profiles,  in  close  agreement  with 
Blasius  profiles,  for  the  quasi-steady  regime.   Very 
recently  King  [13]  has  obtained  similarity  solutions  for 
low  speed  compressible  oscillating  flow  and  for  the 
special  case  of  an  oscillating  wedge  in  a  hypersonic  flow. 

Koob  and  Abbott  [14]  have  recently  developed  a  unique 
integral  method  for  the  analysis  of  the  time  dependent 
laminar  boundary  layer  in  incompressible  flow.   The  method 
is  demonstrated  for  the  case  of  an  accelerating  flat  plate, 
and  shows  promise  of  being  extendable  to  more  complex 
unsteady  flow  problems 

C.   PREVIOUS  EXPERIMENTAL  STUDIES 

In  view  of  the  rather  limited  scope  of  analytic  studies 
of  unsteady  boundary  layer  phenomena,  it  is  somewhat  sur- 
prising that  more  experimental  investigations  have  not  been 
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attempted.   The  relative  lack  of  experimental  data  is 
undoubtedly  due  to  the  problems  encountered  in  establishing 
the  proper  flow  environment  and  accurately  recording  the 
resulting  phenomena. 

Richardson  and  Tyler  [15]  conducted  one  of  the  earliest 
experimental  studies  in  confirming  the  analyses  of  Grace 
[16]  and  Sexl  [17],  who  treated  the  problem  of  oscillating 
flow  in  a  tube  in  the  absence  of  a  mean  flow.   Nickerson 
[7],  in  conjunction  with  his  analytical  work,  made  hot-wire 
measurements  of  the  laminar  boundary  layer  on  a  flat  plate 
oscillating  in  a  Blasius  type  mean  flow  and  obtained  time- 
average  velocity  profiles  and  heat  transfer  data  for  the 
quasi-steady  (low  frequency)  regime.   His  work  confirmed 
in  part,  both  his  and  Ligh thill's  analyses.   Mechanical 
difficulties  precluded  a  more  substantial  experimental 
program.   This  work  was  subsequently  extended  by  Hill  [18] 
who  imposed  oscillations  directly  on  both  Blasius  and 
Howarth  type  mean  flows.   Hill  employed  a  tubular  wind 
tunnel  equipped  with  an  oscillating  gate  valve,  capable 
of  producing  frequencies  of  up  to  four  c.p.s.   In  the  case 
of  the  Howarth  mean  flow,  the  pressure  gradient  was  con- 
trolled by  a  porous  tube  along  the  axis  of  the  test  section, 
Hill's  work  with  the  Blasius  mean  flow  confirmed  both 
Nickerson  and  Lighthill's  analytical  results. 

Hori  [19]  has  investigated  the  laminar  boundary  layer 
on  a  circular  cylinder  oscillating  at  low  frequency  in 
steady  flows  of  extremely  low  velocity.   His  results 


21 


confirm  Lighthill's  quasi-steady  regime  solutions  for 
oscillating  cylinders. 

D.   TRANSITION  AND  SEPARATION 

Relatively  little  work  has  been  reported  on  two  very 
important  aspects  of  unsteady  boundary  layer  flow, 
transition  and  separation.   Many  investigations  have  been 
conducted  on  problems  involving  these  phenomena,  such  as 
"retreating  blade"  stall,  rotating  stall,  flutter,  and 
airfoil  performance  in  the  presence  of  turbulence,  but 
they  have  in  most  cases  been  restricted  to  gross  effects, 
and  rarely  treat  boundary  layer  phenomena. 

Miller  [20]  investigated  transition  phenomena  on  a 
flat  plate  in  a  Blasius  mean  flow  on  which  oscillations 
were  imposed  by  using  a  rotating  shutter  valve.   The 
design  of  the  shutter  valve  allowed  Miller  to  attain 
oscillation  frequencies  from  4  to  125  c.p.s.  and  ampli- 
tudes from  8  to  67  per  cent  of  the  freestream  velocity. 
This  broad  range  of  operating  conditions  makes  this  work 
of  particular  interest  and  significance,  as  prior  experi- 
mental efforts  had  been  mechanically  restricted  to 
frequencies  of  less  than  20  c.p.s.  and  amplitudes  no  greater 
than  10  per  cent  of  the  mean  freestream  velocity.   Miller 
confirms  the  analytical  predictions  of  Liepmann  [21]  and 
Greenspan  and  Benney  [22].   In  addition,  he  has  determined 
the  Transition  Reynolds  Number,  Turbulent  Reynolds  Number 
and  turbulent  intermittency  factor  for  oscillating  Blasius 
type  flow. 
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Even  less  information  exists  on  separation  phenomena 
in  unsteady  flow  than  on  transition.   Analytical  approaches 
of  proven  merit  for  oscillating  Blasius  flows  have  in 
general  proved  to  be  inadequate  for  the  treatment  of 
oscillating  flows  in  the  presence  of  adverse  mean  pressure 
gradients.   Lighthill's  [4]  analysis  of  the  laminar 
boundary  layer  on  a  cylindrical  body  in  oscillating  flow 
loses  its  applicability  soon  after  the  pressure  gradient 
on  the  body  becomes  positive.   Similarly,  Hori ' s  [19] 
analytical  predictions  for  a  cylinder  in  oscillating  flow 
diverge  markedly  from  his  experimental  results  as  soon  as 
adverse  pressure  gradient  effects  are  encountered.   Hill 
[18]  has  analytically  treated  the  case  of  oscillations 
superimposed  on  a  Howarth  mean  flow  and  experimentally 
confirmed  his  solution  for  one  relatively  small  adverse 
pressure  gradient.   He  does  not,  however,  consider  his 
analysis  applicable  in  the  neighborhood  of  separation. 

Both  Hill  and  Hori,  although  conducting  experimental 
investigations  in  adverse  pressure  gradient  regimes, 
present  no  data  relevant  to  laminar  boundary  layer  separ- 
ation.  It  should  be  noted,  however,  that  these  experiments 
were  conducted  in  conjunction  with  analytical  investigations 
and  were  therefore  restricted  to  the  limits  of  applicability 
of  the  analyses.   All  of  the  investigations  mentioned  above 
have  been  primarily  concerned  with  determining  changes 
occurring  in  the  maximum  amplitude  and  phase  relation  of 
the  oscillating  velocity  component  in  the  boundary  layer. 
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It  would  appear  that  knowledge  of  separation  phenomena  can 
most  readily  be  acquired  by  observing  the  changes  occurring 
in  the  complete  velocity  profile  during  a  full  cycle  of 
oscillation. 

Moore  [23,24]  and  more  recently,  Hartunian  and  Moore 
[25] ,  have  made  preliminary  studies  concerned  with  the 
separation  of  the  laminar  boundary  layer  in  unsteady  flow. 
These  studies  indicate  that  the  general  case  of  a  "semi- 
quasi-steady"  oscillating  flow  may  be  divided  into  two 
problems,  first  the  unsteady  boundary  layer  perturbation 
for  a  fixed  separation  point,  and  second,  the  quasi- 
steady  boundary  layer  over  a  surface  which  is  in  motion 
relative  to  a  fixed  separation  point.   These  efforts 
failed  to  attain  their  stated  goal  of  determining  a  mean- 
ingful definition  of  separation  in  unsteady  flow.   Hartunian 
and  Moore  did,  however,  without  rigorous  proof,  predict  the 
general  effects  of  local  accelerations  and  decelerations 
on  the  boundary  layer  profiles,  and  conjectured  that 
separation  might  be  characterized  by  a  simultaneous  vanish- 
of  velocity  and  shear  somewhere  in  the  boundary  layer. 
Perhaps  the  most  significant  feature  of  this  work  is  that 
it  graphically  illustrates  the  difficulty  of  analytically 
arriving  at  a  definition  of  separation  in  unsteady  flow, 
as  well  as  pointing  out  the  need  for  such  a  definition. 
Rott  [26]  also  comments  on  the  problem  of  defining  separ- 
ation in  unsteady  flow  in  a  recent  survey  of  laminar  flow 
phenomena.   Vidal  [27]  attempted  to  experimentally  verify 
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and  to  extend  the  predictions  of  Hartunian  and  Moore. 
Unfortunately,  his  work  was  prematurely  terminated,  and 
although  his  efforts  showed  promise  of  some  success, 
they  must,  on  the  whole,  be  considered  inconclusive. 

The  phenomena  of  laminar  boundary  layer  separation  in 
unsteady  flow  remain  virtually  unexplored,  and  literally 
undefined. 
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II.   LAMINAR  BOUNDARY  LAYER  SEPARATION 
IN  OSCILLATING  FLOW 

A.  INTRODUCTION 

To  some  degree,  all  aerodynamic  flows  are  influenced 
by  unsteadiness.   Commonly  encountered  departures  from 
steady  flow  cover  a  wide  range  of  amplitudes  and  vary 
in  character  from  the  random  perturbations  of  clear  air 
turbulence  to  the  sinusoidal  oscillations  encountered 
by  rotary  wings  and  turbo-machinery  blades.   Flow  phenomena 
involving  periodic  unsteadiness  are  generally  simpler  to 
treat  analytically  and  simulate  experimentally  than 
similar  random  phenomena,  and  any  knowledge  gained  con- 
cerning these  periodic  processes  will  usually  facilitate 
the  study  of  their  random  counterparts.   Consequently, 
the  effects  of  periodic  oscillations  on  aerodynamic  flow 
phenomena  are  of  basic  interest  to  both  researchers  and 
design  engineers.   One  of  the  most  important  of  the 
phenomena  affected  by  unsteadiness  is  laminar  boundary 
layer  separation,  since  its  occurrence  grossly  alters 
the  performance  of  flow  devices.   Thus,  knowledge  of  the 
character  of  laminar  boundary  layer  separation  in  oscillat- 
ing flow  is  of  practical  importance. 

B.  CONCERNING  THE  DEFINITION  OF  SEPARATION 

It  has  been  previously  noted  that  little  is  known 
about  the  behavior  of  the  laminar  boundary  layer  in  an 


26 


oscillating  flow,  adverse  pressure  gradient  environment. 
The  phenomenon  of  boundary  layer  separation  in  such  flows 
remains  virtually  uninvestigated,  and  more  particularly, 
not  well  defined. 

In  steady  flows,  laminar  separation  is  rigorously 
defined  as  commencing  at  that  point  on  the  surface  where 
the  first  indications  of  fluid  convection  in  the  upstream 
direction,  or  "reverse  flow,"  are  observed.   This  reverse 
flow  invariably  forms  eddies  which  decay  and  create  the 
region  of  generally  disturbed  flow  associated  with 
separation.   The  criterion  for  actually  determining  the 
"separation  point,"  consistent  with  the  above  definition, 
is  quite  logically  the  disappearance  of  shear,  or  equiv- 
alently,  velocity  gradient  at  the  wall. 

In  oscillating  flows,  the  steady  flow  description  of 
separation  is  not  applicable.   In  such  flows  the  character 
of  the  boundary  layer  at  any  spatial  point  varies  as  a 
function  of  time  with  the  cyclic  component  of  the  flow. 
Hence,  no  single  velocity  profile  shape  can  be  associated 
with  a  given  point  on  the  surface.   Moore  [2  3]  in  his 
analytical  investigations,  predicts  that  instantaneous 
"zero  wall  shear"  profiles  can  occur  when  the  flow  down- 
stream bears  no  evidence  of  separation.   Stuart  [28]  calls 
attention  to  the  fact  that  in  oscillating  flows,  it  is 
analytically  possible  to  observe  both  instantaneous  "zero 
wall  shear"  and  "reverse  flow"  profiles  even  in  the 
absence  of  any  mean  pressure  gradients.   This  suggests  that 
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regions  of  reverse  flow  occur  periodically  in  oscillating 
boundary  layers  without  the  ensuing  decay  to  disturbed 
flow  as  in  the  steady  flow  case.   It  is  therefore  obvious 
that  the  phenomena  of  laminar  boundary  layer  separation 
must  be  redefined  for  the  case  of  oscillating  flows,  and 
suitable  criteria,  consistent  with  the  redefinition,  must 
be  determined  for  recognizing  the  occurrence  of  separation, 

C.   OBJECTIVES 

The  primary  objective  of  the  present  work  is  to  con- 
duct an  experimental  study  to  determine  a  rational 
definition  of  laminar  boundary  layer  separation,  and 
criteria  for  its  recognition  in  flows  wherein  the  free- 
stream  velocity  is  of  the  form:  U  =  U  +  U  sin  cot. 

The  method  of  attack  consists  of  an  investigation 
into  the  behavior  and  character  of  instantaneous  boundary 
layer  velocity  profiles  in  such  flows  for  as  wide  a 
range  of  freestream  conditions  as  is  practical,  from  which 
a  rational  definition  and  criteria  should  result. 

An  additional  objective  is  to  correlate  the  occurrence 
of  separation  with  flow  parameter  variation  so  that 
reasonable  engineering  predictions  of  separation  in 
oscillating  flows  may  be  made. 
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III.   EXPERIMENTAL  EQUIPMENT 

A.   WIND  TUNNEL 

1.   General  Description 

The  study  was  conducted  in  the  low-speed,  oscillat- 
ing flow  wind  tunnel  located  in  the  Aeronautics  Labora- 
tories of  the  Naval  Postgraduate  School.   This  wind  tunnel 
is  of  open  circuit  design,  with  a  24-inch  square  by  223-inch 
long  test  section.   A  plan  view  of  the  tunnel  is  presented 
in  Figure  1.   The  tunnel  inlet  is  eight  feet  square, 
resulting  in  a  16:1  contraction  ratio.   Three  high  solidity 
screens  located  in  the  inlet  section  just  upstream  of  the 
nozzle  produce  measured  freestream  turbulence  intensities 
of  from  0.261  to  0.413  per  cent  for  the  velocities  encoun- 
tered in  the  present  work. 

The  wind  tunnel  drive  consists  of  two  Joy  Axivane  Fans 
in  series,  each  of  which  has  an  internal,  100  horsepower, 
direct  connected,  1750  rpm  motor.   The  fan  blades  are 
internally  adjustable  through  a  pitch  range  of  25  to  55 
degrees,  providing  a  wide  operating  base.   Two  sets  of 
variable  inlet  vanes,  located  immediately  upstream  of  each 
fan,  are  externally  operated  to  provide  control  of  test 
section  velocity.   These  vanes  are  of  multileaf  design, 
and  preswirl  the  air  in  the  direction  of  fan  rotation  to 
reduce  fan  capacity.   The  total  range  of  tunnel  velocity 
is  from  10  to  250  feet  per  second. 
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2 .   Rotating  Shutter  Valve 

Two  fundamental  methods  of  creating  an  oscillating 
flow  environment  have  been  employed  in  the  past.   Both 
Nickerson  [7]  and  Hori  [19]  introduced  oscillations  by 
oscillating  their  models  in  a  steady  flow  environment. 
This  method  severely  restricts  the  range  of  attainable 
frequencies  because  of  mechanical  complications,  and  also 
introduces  measurement  difficulties.   The  other  approach 
is,  of  course,  to  actually  oscillate  the  flow  over  a 
stationary  model.   Hill  [18]  used  a  sliding  shutter  to 
impose  oscillations  on  the  freestream  but  was  still 
restricted  by  mechanical  difficulties  to  low  frequencies. 

The  most  successful  method  of  obtaining  an  oscillat- 
ing flow  with  large  ranges  of  frequency  and  amplitude  was 
that  employed  by  Karlsson  [29],  and  later  by  Miller  [20] 
in  his  investigation  of  transition.   A  rotating  shutter 
valve,  immediately  downstream  of  the  test  section,  is 
employed  to  superimpose  a  periodic  variation  of  velocity 
on  the  mean  flow.   The  method  used  in  the  present  investi- 
gation is  virtually  identical  to  that  employed  by  Miller. 
The  shutter  valve  consists  of  four  horizontal  steel  shafts 
equally  spaced  across  the  test  section.   The  shafts  are 
slotted  to  accommodate  flat  blades  of  various  widths, 
forming  a  set  of  four  butterfly  valves  spanning  the  test 
section.   Figure  2  is  a  photograph  of  the  shutter  valve. 
Each  blade  is  driven  from  its  immediate  neighbor  by  means 
of  a  timing  belt  and  pulley  arrangement.   The  bottom  shaft 
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is  driven  by  a  five  horsepower  variable  speed  electric 
motor,  through  a  timing  belt  and  pulley.   An  intermediate 
shaft  between  the  motor  and  shutter  valve  permits  a  wide 
variety  of  pulley  ratios.   This  drive  arrangement  provides 
a  frequency  range  of  from  two  cycles  per  second  to  the 
first  critical  frequency  of  9  33  cycles  per  second.   The 
electric  motor  presently  in  use,  however,  restricts  the 
oscillation  frequency  to  a  maximum  of  240  cycles  per  second, 
The  amplitude  of  oscillation  is  controlled  by  blade  width. 
Test  section  closure  may  be  varied  from  25  to  100  per  cent. 
The  resulting  amplitude  of  oscillation  of  test  section 
velocity  is  a  function  of  frequency,  mean  velocity  and 
pressure  gradient.   In  this  investigation,  blades  producing 
50.0,  66.7  and  82.5  per  cent  closure  were  used,  resulting 
in  an  amplitude  range  of  from  3  to  40  per  cent  of  the  local 
mean  freestream  velocity. 
3.   Test  Section 

The  wind  tunnel  test  section  is  shown  in  Figure  3. 
Continuous  pieces  of  two-inch  thick  aluminum,  24  inches 
wide  and  223  inches  long,  form  the  upper  and  lower  test 
section  walls.   Each  of  the  side  walls  consists  of  three 
two-inch  thick  panels  of  stress-relieved  Lucite.   For  this 
investigation,  the  central  side  wall  panel  on  the  opposite 
side  of  the  tunnel  from  the  control  console  was  replaced 
with  two-inch  thick  plywood  to  facilitate  the  mounting  of 
instrumentation.   The  Lucite  panels  on  the  console  side 
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of  the  test  section  are  hinged  and  may  be  raised  hydraul- 
ically,  providing  access  to  the  test  section.   The  heavy 
construction  of  the  test  section  is  intended  to  minimize 
deflections  induced  by  rapid  changes  in  static  pressure. 

Figure  4  is  a  typical  test  section  velocity  pro- 
file. Velocity  variation  is  less  than  one  per  cent  from 
the  mean  to  within  three  inches  of  any  wall. 

Figure  5  is  an  overall  photographic  view  of  the 
wind  tunnel. 

B .   MODEL 

The  model  used  in  this  study  is  an  elongated  two- 
dimensional,  symetric  airfoil  with  a  42  inch  chord  and 
thickness  ratio  of  .0357.   The  model  is  shown  in  Figure  6 
It  was  fabricated  of  mahogany  and,  with  its  supporting 
beams,  spans  the  test  section.   From  the  leading  edge, 
the  first  five  inches  of  the  model  profile  are  a  half 
ellipse  terminating  at  the  minor  axis  which  is  1.5  inches 
long.   The  next  sixteen  inches  consist  of  a  constant 
1.5  inch  thickness  section.   The  last  21  inches  consist 
of  straight  tapers  (smoothly  faired  into  the  preceeding 
section)  from  1.5  inches  to  a  thickness  of  .0625  inches 
on  the  centerline  at  the  trailing  edge.   The  last  twenty 
inches  of  the  model  have  two  three-inch  wide  strips  of 
mild  steel  inlaid  in  the  upper  surface  to  provide  a  chord- 
wise  traversing  surface  for  the  probe  carriage.   These 
strips  are  three  inches  apart,  centered  on  the  mid-span 
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FIGURE     4 

TYPICAL  WIND    TUNNEL  TEST  SECTION    VELOCITY   PROFILE 

24    INCHES  UPSTREAM    OF   MODEL 
UMax  =  20  ft/sec 
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line,  and  extend  1.5  inches  beyond  the  trailing  edge. 
There  are  seventeen  static  pressure  taps  located  in  the 
upper  surface,  distributed  from  leading  to  trailing  edge 
in  a  line  four  inches  from  the  edge  of  the  model. 

Several  features  were  incorporated  in  the  model  design 
to  help  ensure  rigidity  and  prevent  warping.   The  basic 
model  is  of  laminated  construction,  made  from  .75  inch 
thick  slabs  running  spanwise  and  bonded  together  with 
epoxy  over  an  aluminum  spar.   End  plates  made  from  .0625 
inch  thick  aluminum  T-section  were  incorporated,  with 
the  section  base  cut  to  conform  to  the  profile  shape 
and  the  spar  imbedded  in  the  model.   The  relatively 
thin  trailing  edge  was  stiffened  with  a  four-inch  wide, 
.0625  inch  thick  stainless  steel  plate  inserted  spanwise 
along  the  centerline. 

The  model  surface  was  given  multiple  coats  of  lacquer 
and  hand  rubbed  to  a  mirror  finish.   No  discontinuity 
between  the  inlaid  steel  strips  and  the  wooden  surface 
is  detectable. 

For  this  investigation,  the  model  was  mounted  hori- 
zontally on  the  tunnel  centerline  at  zero  angle  of  attack. 

Control  of  the  pressure  distribution  was  accomplished 
by  altering  the  contour  of  the  upper  surface  of  the  test 
section  over  the  model.   This  was  done  by  fastening 
shaped  styrofoam  blocks  to  the  upper  test  section  surface. 
The  contour  of  the  model  surface  assured  the  continuity 
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of  the  first  and  second  derivatives  of  the  chordwise 
pressure  distributions  thus  obtained. 

C.   INSTRUMENTATION 

1.  Freestream  Sensors 

A  conventional  pitot-static  tube  and  a  fixed 
hot-wire  anemometer  probe  were  located  twelve  inches 
forward  of  the  model  leading  edge  five  inches  above  the 
centerline  and  were  used  to  determine  mean  freestream 
velocity.   Dynamic  pressure  was  read  from  a  Meriam 
micro-manometer.   The  hot  wire  probe  was  used  to  record 
freestream  turbulence  level  and  to  make  fine  adjustments 
to  the  mean  velocity  settings. 

The  frequency  of  oscillation  was  obtained  from 
a  magnetic  pickup  mounted  outboard  of  the  uppermost 
shutter  blade  shaft,  as  seen  in  Figure  2.   The  output 
frequency  was  read  on  a  Berkeley  decade  counter. 

2 .  Hot-Wire  Anemometers 

Two  constant  temperature  hot-wire  anemometer 
units  were  used  in  this  investigation.   The  primary  unit 
was  the  ten  channel  Security  Associates  Model  1000,  shown 
in  Figure  7.   A  single  channel  Security  Associates  Model 
100  hot-wire  anemometer  was  also  used  as  a  "mobile"  unit. 
It  was  thus  possible  to  operate  and  monitor  eleven  inde- 
pendent hot-wire  sensors  simultaneously.   All  eleven 
basic  circuits,  exclusive  of  power  supplies,  potentiometers 
and  switches,  are  identical  and  interchangeable  solid 
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state  printed  circuit  plug-in  modules.   The  circuits  in 
the  ten  channel  unit  have  common  power  supplies  and 
volt- ammeter .   An  important  feature  of  these  units  is 
that  each  hot-wire  bridge  circuit  is  coupled  to  an  ana- 
logue computer  circuit  which  inverts  King's  Law  and  pro- 
vides an  output  which  is  linear  in  velocity  and  zero  for 
zero  velocity. 

The  "wires"  used  as  sensing  elements  were  0.00015 
inch  diameter  tungsten  filaments,  .125  inches  long.   The 
ends  of  the  filaments  were  copper  plated  to  facilitate 
mounting  on  the  probes  (described  in  the  next  section) 
resulting  in  an  effective  sensing  length  of  from  .080  to 
.090  inches.   The  plating  procedure  and  equipment  is 
described  in  Reference  [30]. 

The  linearity  of  each  of  the  circuits  used  was 
checked  in  steady  freestream  conditions  using  the  pitot- 
static  tube  and  freestream  hot-wire  probe  referred  to 
above.   Figure  8  is  a  typical  calibration  curve.   Response 
at  typical  boundary  layer  velocities  was  further  checked 
by  obtaining  laminar  boundary  layer  velocity  profiles  on 
the  flat  plate  model  described  in  Reference  [31]  and 
comparing  the  results  to  the  Blasius  profile.   Figure  9 
is  typical  of  these  results. 

Since  the  hot-wire  circuits  are  D.C.  coupled, 
total  instantaneous  velocity  was  available  at  the  out- 
puts.  The  direct  current  components  of  the  outputs, 
which  are  proportional  to  the  mean  velocity  components, 
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TYPICAL  HOT  WIRE  ANEMOMETER  CALIBRATION  CURVE 
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TYPICAL     HOTWIRE  ANEMOMETER    CALIBRATION 
VELOCITY    PROFILE   IN  BLASIUS    FLOW 
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were  measured  with  the  precision  voltmeters  contained  in 
the  two  hot-wire  anemometer  units.   The  alternating  current 
components,  proportional  to  the  oscillating  velocity  com- 
ponents, were  displayed  on  a  Tektronix  555  Dual-Beam 
oscilloscope  equipped  with  a  pair  of  4-channel  preamplifier 
units.   Thus,  of  the  eleven  channels  being  monitored, 
eight  of  the  oscillating  velocity  component  outputs  could 
be  simultaneously  displayed.   The  composition  of  the  display 
could  be  easily  changed  by  using  a  switching  unit  incor- 
porated in  the  output  circuitry.   The  resulting  displays 
were  photographically  recorded  with  a  Tektronix  oscillo- 
scope camera. 

The  output  signal  from  the  shutter  blade-actuated 
magnetic  pickup  was  used  as  an  external  trigger  to  control 
oscilloscope  sweep  timing. 

Figure  10  is  a  schematic  diagram  of  the  instru- 
mentation circuitry.   Figure  11  is  a  photograph  of  the 
control  console,  showing  the  basic  instrumentation  units. 

D.   MULTIPLE  HOT-WIRE  PROBE 

1 .   General  Characteristics 

The  prime  requirement  for  the  hot-wire  probe  in 
this  investigation  was  to  place  a  chordwise-traversable 
array  of  eleven  hot-wires  in  the  boundary  layer.   Due  to 
the  oscillating  character  of  the  flow  environment,  the 
elimination  of  vibration-induced  relative  motion  between 
the  hot-wires  and  the  model  surface  was  an  important 
consideration . 
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2  .   Probe  Modules 

Initial  attempts  at  probe  design  centered  on  the 
concept  of  having  ten  hot-wire  filaments,  parallel  to 
each  other  and  the  model  surface,  oriented  in  a  single 
vertical  array  of  a  depth  corresponding  to  expected 
boundary  layer  thicknesses  (approximately  .25  to  .35  inches). 
Such  an  array  proved  to  be  physically  impractical  to  con- 
struct.  Even  when  the  number  of  wires  in  the  array  was 
reduced  to  five,  the  mounting  of  the  hot-wire  filaments  to 
the  probe  needles  proved  to  be  extremely  difficult  due  to 
the  proximity  of  the  wires.   In  addition,  these  arrays 
required  that  the  probe  needles  used  be  of  such  small 
diameter,  and  so  flexible,  that  they  would  almost  cer- 
tainly vibrate  in  response  to  the  oscillating  velocities. 

It  was  finally  decided  to  divide  the  ten  wire 
array  into  five  two-wire  subarrays  extending  in  a  spanwise 
direction.   Each  two-wire  subarray  is  mounted  on  a  separate 
plug-in  probe  module,  as  shown  in  Figure  12.   Each  module 
consists  of  a  small  block  of  epoxy  resin,  in  which  the 
four  necessary  probe  needles  and  their  associated  male 
electrical  contacts  are  cast.   Jeweler's  broaches  were 
used  for  probe  needles.   They  extend  at  a  45  degree  angle 
relative  to  the  module  block.   The  tips  are  spaced  .100 
inches  apart  vertically  and  .150  inches  apart  horizontally 
in  a  plane  normal  to  the  model  surface  when  the  module  is 
mounted  in  the  probe  carriage.   The  hot-wire  filaments  were 
soldered  to  the  probe  tips,  with  the  aid  of  a  stereo 
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microscope  and  specially  designed  holding  jig,  so  that  the 
filaments  are  mutually  parallel.   In  order  to  eliminate 
possible  displacement  of  the  probe  needles  relative  to 
each  other,  an  epoxy  web  was  formed  between  each  two 
vertically  adjacent  needles.   Ten  of  these  probe  modules 
were  constructed  and  used  in  the  investigation. 

3 .  Vertical  Traverse  Probe 

A  vertically  traversable  probe,  controlled  by  an 
external  micrometer  through  a  flexible  shaft  arrangement, 
was  also  incorporated  in  the  hot-wire  array.   Figure  13 
shows  the  traverse  probe  disassembled,  and  Figure  14,  the 
micrometer  control. 

4 .  Probe  Carriage 

The  fixed  array  probe  modules  and  the  vertical 
traverse  probe  were  both  mounted  on  a  chordwise-traversable 
aluminum  probe  carriage.   Figure  15  is  a  photograph  of 
the  carriage  with  probes  mounted.   The  carriage  consists 
basically  of  an  aerodynamically  shaped  crosspiece  housing 
the  probes  and  vertical  support  legs.   The  crosspiece  is 
provided  with  twelve  spanwise  plug-in  stations  for  the 
modules,  each  station  having  a  vertical  set  screw  which 
is  used  to  vary  the  position  of  the  hot-wires  in  the  array 
plane.   This  arrangement  provides  a  maximum  fixed  array 
depth  of  0.40  inches  normal  to  the  model  surface.   The 
spanwise  width  of  the  array  may  be  varied  from  1.4  to  4.6 
inches  by  using  various  combinations  of  the  plug-in  sta- 
tions.  The  vertical  traverse  probe  assembly  is  mounted 


50 


'• 


m,    « 


■""'?. 


.. 


Q 
UJ 

_J 
CD 

UJ 

in 

o 


ro 

LU 

a: 

3 
O 
Li. 


UJ 
GO 
O 

q: 
q_ 

in 

LU 


< 
O 

UJ 
> 


UJ 

X 


X 

a. 
< 
a: 

o 

& 

x 


51 


Ll! 

<r 


_J 
O 


O 
O 

a: 

Ll! 
H 
LlI 

o 

LlI 
(/) 
(T. 
LU 

1 

_l 
< 
O 

H 

a: 

UJ 

X 

h- 

li. 
o 

X 
CL 
< 

a: 

g 

o 

X 
Q_ 


52 


FIGURE    15     , 

PHOTOGRAPH    OF  THE    PROBE    CARRIAGE 
WITH    MODULES    AND  VERTICAL   TRAVERSE  PROBE 

MOUNTED 
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on  the  crosspiece  center  line  such  that  its  hot-wire  always 
lies  in  the  same  plane  normal  to  the  model  surface  as  the 
fixed  array  wires. 

Three  teflon-coated  Alnico  sintered  magnets  are 
mounted  to  the  outboard  side  of  each  support  leg  at  its 
base.   The  magnets  hold  the  probe  carriage  in  positive 
contact  with  the  steel  inlays  in  the  model,  while  still 
allowing  the  probe  assembly  to  traverse  in  the  chord- 
wise  direction. 

The  structure  of  the  probe  carriage,  coupled 
with  the  magnetic  mode  of  attachment  to  the  model,  virtually 
eliminates  the  possibility  of  oscillation  induced  relative 
motion  between  the  probes  and  the  model  surface. 

The  possibility  of  boundary  layer  profile  distortion 
in  the  array  plane  due  to  the  proximity  of  the  probe 
carriage  was  investigated  in  a  steady  flow,  positive 
pressure  gradient  environment.   Figure  16  shows  the  results 
of  this  investigation.   Apparent  velocity  distortions  of 
up  to  five  per  cent  were  observed.   Since  the  distortions 
were  consistent  across  the  array,  and  had  little  effect 
on  the  general  character  of  the  profile,  they  were  con- 
sidered acceptable  as  a  necessary  part  of  the  experimental 
environment. 

The  probe  carriage  was  moved  fore  and  aft  on  the 
model  by  means  of  a  remotely  controlled  traversing  mech- 
anism mounted  on  the  test  section  centerline  aft  of  the 
model.   The  chordwise  traverse  range  covers  the  aftermost 
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FIGURE    16 
EFFECT  OF  PRESENCE  OF  PROBE  CARRIAGE  ON  VELOCITY  PROFILE 
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21  inches  of  the  model.   A  view  of  the  chordwise  traverse 
mechanism  in  place  is  shown  in  Figure  17.   Figure  18  shows 
the  probe  assembly  and  model  mounted  in  the  test  section. 
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IV.   EXPERIMENTAL  PROCEDURE 

A.   CHARACTER  OF  THE  FLOW  ENVIRONMENT 

Prior  to  the  collection  of  experimental  data,  it  was 
necessary  to  investigate  the  nature  of  the  flow  field. 
The  use  of  the  spanwise  distributed  hot-wire  array  to 
observe  boundary  layer  phenomena  assumes  the  flow  in  the 
regime  of  interest  to  be  two-dimensional.   To  determine 
if  this  was  indeed  the  case,  all  of  the  array  probes  were 
set  at  the  same  level  above  the  model  surface.   The  probe 
assembly  was  then  traversed  through  its  full  chordwise 
range  for  various  oscillating  flow  conditions.   This  pro- 
cedure was  repeated  for  various  vertical  displacements 
spanning  the  boundary  layer  and  extending  into  the 
adjacent  freestream.   Comparison  of  output  displays  for 
the  several  wires  verified  that  the  flow  was  indeed  two- 
dimensional  throughout  the  experimental  regime. 

Further  preliminary  surveys,  made  with  the  arrays  of 
probes  set  at  varying  depths,  revealed  that  the  flow  for 
all  conditions  of  interest  was  "strictly  periodic. "   As 
used  here,  "strictly  periodic"  refers  to  the  fact  that  at 
any  spatially  fixed  position  in  the  flow,  the  wave  form 
representing  the  oscillating  component  of  velocity 
exhibited  no  aperiodic  changes  in  its  amplitude,  phase, 
or  shape. 

The  strictly  periodic  quality  of  the  flow  made  possible 
the  employment  of  relatively  simple,  and  consequently 
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accurate,  methods  in  collecting  data.   While  the  fixed 
probe  array  provided  data  that  could  be  used  to  obtain 
truly  instantaneous  boundary  layer  information,  the 
relatively  small  number  of  data  points  available  (8) , 
and  the  uncertainty  in  determining  probe  depth  in  the 
region  immediately  adjacent  to  the  wall,  restricts  the 
detail  and  accuracy  of  the  profiles  obtained.   However, 
since  the  flow  was  strictly  periodic,  the  vertical 
traverse  probe  could  be  employed  to  provide  very  accurate 
profiles,  that,  while  not  in  themselves  instantaneous, 
are  in  every  way  equivalent  to  true  instantaneous  profiles. 
The  fixed  array  was  still  of  great  value  in  qualitatively 
observing  the  instantaneous  behavior  of  the  boundary  layer, 
as  well  as  in  supplementing  and  verifying  observations  made 
with  the  vertical  traverse  probe. 

B.   INVESTIGATIVE  PROCEDURE 

Preliminary  investigations  were  carried  out  by 
traversing  the  entire  adverse  pressure  gradient  regime 
with  the  probe  assembly.   The  output  of  the  fixed  probe 
array  was  monitored  on  the  anemometer  voltmeter  and  the 
oscilloscope  display,  providing  a  continuous  qualitative 
picture  of  boundary  layer  behavior  for  a  wide  variety  of 
flow  conditions.   When  regions  of  particular  interest 
were  observed,  accurate  "instantaneous"  boundary  layer 
profiles  were  obtained  by  using  the  vertical  traverse 
probe.   For  each  point  in  the  traverse,  the  steady  (mean) 
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velocity  component  of  the  flow  was  recorded  from  the  ane- 
mometer voltmeter,  and  the  oscillating  component  was 
recorded  photographically  from  the  oscilloscope  display. 
For  these  photographs,  the  display  was  composed  of  two 
sweeps,  the  uppermost  presenting  the  output  of  an  array 
hot-wire  positioned  in  the  freestream,  and  the  lower  the 
output  of  the  traverse  probe.   Boundary  layer  profiles 
could  then  be  plotted  by  vectorially  adding  the  mean 
velocities  for  each  point  in  the  traverse  with  the  appro- 
priate oscillating  component,  measured  from  the  photographs 
at  the  desired  point  in  the  cycle. 

C.   FLOW  VARIABLES 

1.  Mean  Velocity 

The  initial  profiles  acquired  indicated  that 
boundary  layer  transition  effects  were  encountered  for 
mean  flow  velocities  of  30  feet  per  second  and  greater. 
At  mean  velocities  less  than  15  feet  per  second,  the  pro- 
files were  highly  susceptable  to  disturbances  caused  by 
atmospheric  wind  gusts  at  the  tunnel  entrance.   Therefore, 
mean  velocities  of  15,  20  and  25  feet  per  second  were  used 
in  the  investigation. 

2 .  Static  Pressure  Distribution 

Three  different  mean  static  pressure  distributions 
were  employed  in  the  investigation.   Figure  19  is  a  chord- 
wise  plot  of  the  resulting  mean  pressure  coefficient  dis- 
tributions . 
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MEAN     PRESSURE     COEFFICIENT   DISTRIBUTIONS 
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All  three  distributions  used  were  designed  to  pro- 
vide a  significant  chordwise  region  of  favorable  pressure 
gradient  upstream  of  the  beginning  of  the  adverse  pressure 
gradient  region.   This  was  done  to  introduce  separation 
inhibiting  qualities  into  the  boundary  layer,  so  that  the 
separation  process  would  be  spatially  extended,  facilitat- 
ing observation  of  the  relevent  profile  phenomena.   The 
magnitude  of  the  favorable  pressure  gradient,  and  its 
proximity  to  the  adverse  pressure  gradient  regime  were 
varied  in  each  of  the  three  cases  to  provide  a  wide  range 
of  the  pressure  parameter.   Values  of  the  dimensionless 
pressure  parameter  (defined  below)  were  graphically  obtained 
from  the  dimensional  forms  of  the  distributions  in  Figure 
19. 

3 .   Frequency  of  Oscillation 

Early  in  the  investigation  it  became  evident  that 
for  the  mean  velocities  used,  three  attainable  frequencies 
produced  wave  forms  that  came  as  close  as  possible  to 
being  sinusoidal,  and  were  basically  symmetric.   These 
were  74,  142,  and  226  cycles  per  second. 

Harmonic  analyses  conducted  by  Laferty  and  Hicks 
[33]  on  similar  wave  forms,  produced  by  this  equipment, 
indicate  that  the  average  contributions  of  second  and 
third  harmonics  are  respectively  2  3  and  4  per  cent  of  the 
fundamental.   Similar  analyses  performed  by  Miller  [30] 
indicate  that  four  terms  of  the  Fourier  series  adequately 
represent  this  type  waveform  to  within  one  per  cent. 
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4 .   Amplitude  of  Oscillation 

Shutter  blades  resulting  in  87.5,  66.7,  and  50.0 
per  cent  tunnel  closure  were  used,  producing  amplitudes 
of  oscillation  which  varied  from  40  to  3  per  cent  mean 
freestream  velocity. 

D.   DIMENSIONLESS  FLOW  PARAMETERS 

Commonly  used  forms  of  dimensionless  velocity,  fre- 
quency and  amplitude  were  employed  in  this  investigation, 
and  are  defined  in  the  List  of  Symbols.   The  dimensionless 
pressure  parameter  used,  is  however,  to  the  best  of  this 
author's  knowledge,  unique  to  this  work. 

The  following  parameter  was  used  to  characterize  the 

mean  pressure  gradient  effects  on  the  boundary  layer: 
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The  first  term  in  the  product  accounts  for  the  relative 
effects  of  the  spatial  history  of  the  mean  pressure 
gradient  distribution.   It  is  known  that  the  presence  of 
a  regime  of  favorable  pressure  gradient  upstream  of  the 
adverse  pressure  gradient  tends  to  delay  separation.   The 
extent  of  this  delay  is  increased  as  the  magnitude  of  the 
favorable  pressure  gradient  and  its  proximity  to  the 
adverse  pressure  gradient  regime  increases.   With  this  in 
mind,  the  first  term  was  formulated  as  the  percentage  of 
overall,  spatially  weighted,  pressure  gradient  distribution 
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which  is  adverse,  or  separation  inducing.   For  the  special 
case  of  Howarth-type  flow,  this  history  term  becomes  equal 
to  unity,  due  to  the  absence  of  a  favorable  pressure 
gradient  regime. 

The  complete  pressure  parameter  accounts  therefore, 
for  both  the  magnitude  of  the  local  pressure  gradient 
at  separation  and  the  influence  of  the  preceeding  pressure 
gradient  distribution. 
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V.   RESULTS 

A.   FLOW  PHENOMENA 

1 .   Flow  Region  Classification 

Preliminary  surveys  made  with  the  multiple  probe 
array  revealed  that  the  flow  field  produced  by  the  adverse 
pressure  gradient  could  be  divided  into  three  distinct 
regions.   Proceeding  in  a  streamwise  direction  the  first 
of  these  regions  was  characterized  by  the  regularity  and 
repetitive  nature  of  the  oscillating  velocity  waveform 
in  the  boundary  layer.   In  this  region,  random  pertur- 
bations in  the  waveform  were  infrequent,  and  when  observed, 
were  quickly  damped.   When  velocity  fluctuations  of  this 
nature  occurred,  it  was  only  in  conjunction  with  similar 
perturbations  in  the  freestream  velocity.   This  portion  of 
the  flow  may  be  termed  the  "pre-wake  region."   The  pre- 
wake  region  ends  at  a  downstream  point  where  random  per- 
turbations on  the  oscillating  velocity  waveform  begin  to 
originate  in  the  boundary  layer,  persist,  and  grow  with 
downstream  displacement.   No  associated  velocity  fluc- 
tuations are  observed  in  the  freestream  flow.   This  region 
will  be  referred  to  as  the  "region  of  wake  formation." 
Velocity  perturbations  in  the  region  of  wake  formation 
increase  in  severity  until  they  finally  begin  to  induce 
significant  random  fluctuations  in  the  local  freestream 
velocity,  marking  the  beginning  of  the  "region  of  wake 
flow. "   Figure  20  presents  a  sequence  of  waveform 
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FIGURE    20 

TYPICAL  VELOCITY  WAVEFORMS  IN  THE  THREE  ADVERSE 
PRESSURE  GRADIENT  FLOW  REGIONS 
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photographs  along  a  boundary  layer  mean  streamline,  illus- 
trating the  characteristics  of  these  three  distinct  flow 
regions.   In  these  photographs,  as  well  as  subsequent 
waveform  sequences,  the  upper  oscilloscope  sweep  depicts 
the  oscillating  velocity  component  in  the  local  freestream, 

Also,  in  all  waveforms  presented,  the  voltage,  or 
velocity  axis  is  reverse-oriented,  such  that  velocity 
increases  progress  in  the  "downward"  or  "minus-y"  co- 
ordinate direction.   This  is  due  to  voltage  inversion 
inherent  in  the  anemometer  circuit.   The  time  coordinate 
is  oriented  in  the  usual  manner,  i.e.,  increasing  from 
left  to  right. 

2 .   Flow  Reversal 

Multiple  probe  array  observations  revealed  that 
certain  periodically  occurring  phenomena  were  consistently 
present  in  the  boundary  layer  velocity  waveforms.   The 
most  significant  of  these  effects  was  characterized  by  an 
initial  flattening  of  the  waveform  peak  at  the  position  of 
minimum  local  velocity.   This  first  occurs  in  the  lower 
half  of  the  boundary  layer.   As  the  hot-wire  is  moved 
closer  to  the  wall,  the  distortion  intensifies,  resulting 
in  an  apparent  dip  in  the  waveform,  indicating  a  transient 
relative  increase  in  velocity.   This  phenomenon  was  inter- 
preted as  an  indication  of  the  existence  of  reverse  flow 
during  the  affected  portion  of  the  cycle.   This  interpre- 
tation was  based  on  the  fact  that  the  hot-wire  anemometer 
is  an  absolute  value  sensor,  hence  for  an  instantaneous 


68 


reverse  flow  condition,  the  anemometer  will  present  the 
negative,  or  reversed,  velocity  portion  of  the  waveform  as 
positive,  resulting  in  a  "folded"  waveform.   This  deduction 
was  unquestionably  verified  by  plotting  the  instantaneous 
profiles  for  that  portion  of  the  cycle,  and  observing  that 
the  waveform  dip  must  be  considered  as  negative,  or  reversed, 
velocity  to  preserve  the  continuity  of  the  velocity  deriv- 
atives.  Figure  21  presents  a  typical  waveform  sequence, 
observed  during  a  vertical  traverse  of  the  boundary  layer, 
in  which  the  reversal  phenomena  is  evident. 

Prior  to  measuring  the  oscillatory  velocity  compo- 
nents as  presented  in  the  oscilloscope  photographs,  the 
areas  enclosed  by  the  waveform  and  the  calibrated  reference 
line  were  graphically  integrated.   From  this  integration, 
corrections  to  the  mean  velocity  reference  position  on 
the  oscilloscope  grid  were  determined  and  applied  to  the 
oscillatory  velocity  component  measurements.   The  oscilla- 
tory components  were  measured  relative  to  this  corrected 
mean  velocity  reference  as  long  as  no  reversal-dip  was 
evident  in  the  waveform.   Since  the  hot-wire  anemometer 
is  an  absolute  value  sensor,  the  occurrence  of  flow 
^reversal  introduces  an  additive  error  in  the  mean  velocity 
read  from  the  anemometer  voltmeter.   Hence,  for  waveforms 
having  a  reversal-dip,  the  nodes  on  either  side  of  the 
dip  were  recognized  as  establishing  a  reference  coincident 
with  zero  total-instantaneous  velocity,  and  measurements 
were  then  made  relative  to  this  reference,  recognizing 
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that  the  "folded"  portion  of  the  waveform  represented 
negative  velocity. 

Considering  anemometer  accuracy  and  the  measurement 
procedure,  values  of  instantaneous  velocity  obtained  in 
the  above  manner  are  estimated  to  have  an  uncertainty 
of  five  per  cent.   Due  to  uncertainties  in  the  determin- 
ation of  mean  boundary  layer  thickness,  the  dimensionless 
normal  spatial  coordinate  values  have  an  estimated 
uncertainty  of  two  per  cent. 
3 .   Periodic  Effects 

In  addition  to  the  reversal-dip,  two  other  period- 
ically occurring  phenomena  were  observed  in  the  boundary 
layer  velocity  waveforms.   These  were  evident,  in  varying 
degrees,  for  most  operating  conditons,  and  are  charac- 
terized as  follows: 

"Amplif ication"-the  growth,  or  amplification  of 
the  oscillatory  waveform  beyond  its  local  free- 
stream  magnitude  as  the  boundary  layer  is 
vertically  traversed; 

"Phase  Shift"-a  phase  lead,  or  tendency  of  the 
oscillatory  changes  in  the  boundary  layer  to 
occur  earlier  in  time  than  respective  changes  in 
the  freestream  flow. 
Both  of  these  effects  may  be  observed  in  Figure  21.   Amp- 
lification in  this  case  is  quite  pronounced,  and  a  maximum 
phase  lead  of  40  degrees  is  observed  near  the  surface. 
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The  magnitude  and  frequency  of  occurrence  of  these 
effects  will  be  discussed  in  a  subsequent  section.   They 
are  introduced  at  this  point  because  of  their  appearance 
in  some  of  the  results  presented. 

B.   INSTANTANEOUS  VELOCITY  PROFILES 

A  total  of  117  data  runs  were  made  to  obtain  instan- 
taneous velocity  profile  information.   Of  these,  the  first 
35  were  made  for  purposes  of  familiarization,  and  to 
determine  practical  limits  of  the  operating  parameters. 
The  remaining  82  runs  yielded  velocity  profile  sequences 
used  in  the  investigation  of  separation  phenomena.   The 
operating  parameters  for  these  latter  runs  are  listed 
in  Table  I,  and  representative  numerical  data  for  each  run 
are  presented  in  Appendix  A. 

Figure  22  illustrates  typical  behavior  of  an  instan- 
taneous velocity  profile  at  a  point  in  the  pre-wake  region, 
during  a  half  cycle  of  oscillation.   Figure  23  presents 
the  waveform  photographs  from  which  the  instantaneous 
velocity  profiles  were  obtained.   It  should  be  noted  that 
in  Figure  22,  as  well  as  in  subsequent  profile  sequences 
to  be  presented,  the  instantaneous  boundary  layer  velocity 
is  non-dimensionalized  with  respect  to  the  local  mean  free- 
stream  velocity,  so  that  temporal  profile  variations  are 
readily  apparent. 

The  profile  sequence  in  Figure  22  was  obtained  well 
upstream  of  the  region  of  wake  formation.   It  is  apparent 
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TABLE  I 

OPERATING  CONDITIONS  FOR  INSTANTANEOUS 
VELOCITY  PROFILE  RUNS 


RUN 

N_  xlO'5 
Re 

Nfxl05 

NA 

x/c 

1 

2.745 

2.893 

.32 

.612 

2 

2.745 

5.552 

.28 

.612 

3 

2.745 

8.837 

.19 

.612 

4 

2.775 

2.893 

.25 

.619 

5 

2.775 

5.552 

.30 

.619 

6 

2.775 

8.837 

.25 

.619 

7 

2.809 

2.893 

.28 

.626 

8 

2.809 

5.552 

.31 

.626 

9 

2.809 

8.837 

.27 

.626 

10 

2.831 

2.893 

.31 

.631 

11 

2.831 

5.552 

.31 

.631 

12 

2.831 

8.837 

.26 

.631 

13 

2.841 

2.893 

.30 

.633 

14 

2.841 

5.552 

.35 

.633 

15 

2.841 

8.837 

.30 

.633 

16 

2.155 

5.144 

.29 

.643 

17 

2.292 

5.144 

.22 

.662 

18 

2.292 

9.871 

.32 

.662 

19 

2.292 

15.709 

.22 

.662 

20 

2.767 

2.893 

.28 

.619 
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TABLE    I     (continued) 


RUN 

NRexlO"5 

Nfxl05 

NA 

x/c 

21 

2.745 

2.893 

.29 

.612 

22 

2.703 

2.893 

.32 

.605 

23 

2.671 

2.893 

.32 

.598 

24 

2.639 

2.893 

.35 

.590 

25 

2.745 

5.552 

.35 

.612 

26 

2.767 

5.552 

.37 

.619 

27 

2.799 

5.552 

.35 

.626 

28 

2.831 

5.552 

.34 

.633 

29 

2.863 

5.552 

.33 

.640 

30 

2.895 

5.552 

.35 

.648 

'  31 

3.023 

8.837 

.24 

.676 

32 

2.991 

8.837 

.34 

.669 

33 

2.959 

8.837 

.36 

.662 

34 

2.927 

8.837 

.39 

.655 

35 

2.903 

8.837 

.36 

.648 

36 

2.863 

8.837 

.38 

.640 

37 

2.863 

5.552 

.36 

.640 

38 

2.852 

5.552 

.34 

.636 

39 

2.852 

5.552 

.34 

.636 

40 

2.809 

5.552 

.36 

.626 

41 

2.075 

9.871 

.35 

.619 

42 

2.099 

9.871 

.35 

.626 

43 

2.131 

9.871 

.37 

.636 

44 

2.155 

9.871 

.35 

.643 

45 

2.203 

9.871 

.34 

.657 
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TABLE  I  (continued) 


RUN 

M_  xlO  5 
Re 

Nfxl05 

NA 

x/c 

46 

2.235 

9.871 

.34 

.667 

47 

1.996 

5.144 

.32 

.595 

48 

2.020 

5.144 

.35 

.602 

49 

2.052 

5.144 

.36 

.612 

50 

2.076 

5.144 

.37 

.619 

51 

2.108 

5.144 

.34 

.629 

52 

2.148 

5.144 

.32 

.640 

53 

2.474 

15.709 

.25 

.738 

54 

2.451 

15.709 

.25 

.731 

55 

3.373 

3.553 

.28 

.610 

56 

3.341 

3.553 

.36 

.602 

57 

3.259 

1.853 

.30 

.583 

58 

3.309 

1.853 

.32 

.595 

59 

3.459 

5.655 

.30 

.619 

60 

3.373 

5.655 

.34 

.610 

61 

2.155 

5.144 

.05 

.643 

62 

2.195 

5.144 

.06 

.655 

63 

2.179 

5.144 

.06 

.650 

64 

2.219 

5.144 

.04 

.667 

65 

2.235 

5.144 

.04 

.667 

66 

2.283 

9.871 

.04 

.681 

67 

2.321 

15.709 

.03 

.690 

68 

2.826 

2.893 

.05 

.631 

69 

2.852 

5.552 

.03 

.638 

70 

2.890 

8.837 

.03 

.645 
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TABLE  I  (continued) 


RUN 

ND  xlO'5 
Re 

Nfxl05 

NA 

x/c 

71 

3.405 

1.853 

.04 

.607 

72 

3.421 

3.553 

.04 

.612 

73 

3.453 

5.655 

.03 

.617 

74 

2.187 

5.144 

.10 

.650 

75 

2.254 

9.871 

.10 

.671 

76 

2.297 

15.709 

.08 

.686 

77 

2.801 

2.893 

.14 

.624 

78 

2.864 

5.552 

.09 

.636 

79 

2.903 

8.837 

.07 

.640 

80 

3.309 

1.853 

.09 

.593 

81 

3.325 

3.553 

.08 

.598 

82 

3.357 

5.655 

.09 

.600 

All  runs  in  Table  I  were  conducted  in  Pressure  Distribution  I. 
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FIGURE  22 

TYPICAL   SEQUENCE  OF  INSTANTANEOUS 
PROFILES  DURING  HALF-CXLE     OF 
OSCILLATION  .    PRE-WAKE  REGION . 

WEAK   TRANSIENT  FLOW  REVERSAL 
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FIGURE    23 


OSCILLATING     VELOCITY     COMPONENT    WAVEFORMS 
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that  a  reverse  flow  condition  exists  during  approximately 
one-half  of  the  full  cycle  of  oscillation,  reaching  its 
maximum  reversal  when  the  instantaneous  freestream  velocity 
is  at  its  minimum  in  the  cycle.   Transient  flow  reversal  of 
this  nature  was  found  to  be  present  to  some  extent  in  all 
adverse  pressure  gradients.   For  amplitudes  of  oscillation 
of  less  than  eight  per  cent  of  the  local  mean  freestream 
velocity,  transient  flow  reversal  first  occurs  at  some 
point  in  the  pre-wake  region.   For  larger  amplitudes, 
transient  flow  reversal  was  observed  throughout  the  pre- 
wake  region.   This  last  behavior  suggests  that  transient 
flow  reversal  may  even  occur  in  the  zero  pressure  gradient 
region,  as  predicted  by  Stuart  [28]  and  Moore  [23]  for 
sufficiently  large  amplitudes.   However,  being  beyond  the 
scope  of  the  present  work,  this  was  not  investigated. 

The  profile  sequence  in  Figure  22  illustrates  a  case 
of  relatively  weak  transient  flow  reversal.   For  purposes 
of  comparison,  Figure  24  illustrates  the  nature  of  a  strong 
transient  flow  reversal.   Figure  25  presents  the  associated 
waveforms.   It  can  be  seen  in  Figures  22  and  24  that  the 
maximum  and  minimum  velocity  profiles  form  an  envelope 
within  which  all  profile  variation  is  contained. 

Moving  downstream,  profile  inflection  gradually  became 
more  pronounced  throughout  the  cycle.   The  duration  of 
transient  flow  reversal  in  the  cycle,  however,  increased 
only  slightly  with  downstream  displacement.   At  a  point 
generally  between  two  and  four  boundary  layer  thicknesses 
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FIGURE  24 

TYPICAL    SEQUENCE  OF  INSTANTANEOUS 
PROFILES  DURING    HALF-CYCLE   OF 
OSCILLATION.  PRE-WAKE  REGION. 

STRONG  TRANSIENT  FLOW  REVERSAL 
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upstream  of  the  region  of  wake  formation,  the  reversed 
flow  portion  of  the  cycle  began  to  expand  rapidly  as  the 
"least  reversed"  profile  approached  a  zero  velocity 
gradient  condition  at  the  wall.   Finally  at  a  point  within 
half  a  boundary  layer  thickness  of  the  beginning  of  the 
region  of  wake  formation,  the  entire  cycle  of  profiles 
exhibited  reversed  flow  or  zero  velocity  gradient  at  the 
wall.   This  behavior  was  common  to  all  operating  conditions, 
and  varied  only  in  the  magnitude  of  the  initial  cyclic 
fraction  of  reversed  flow.   Figure  26  illustrates  the  growth 
of  the  reverse  flow  portion  of  the  cycle  for  several 
operating  conditions.   Figure  27  illustrates  typical  profile 
envelope  changes  with  increasing  transient  flow  reversal. 

Some  mention  should  be  made  of  the  influence  of 
amplification  effects  and  phase  leads  on  the  instantaneous 
profile  sequences.   Amplification  merely  tends  to  locally 
broaden  the  profile  envelope,  without  affecting  the 
general  profile  shape,  or  the  occurrence  of  transient  flow 
reversal.   Conversely,  phase  leads  tend  to  cause  a  local 
narrowing  of  the  profile  envelope,  in  addition  to  that 
caused  by  reducing  freestream  amplitudes.   They  have  little 
influence  on  general  profile  shape,  even  large  phase  leads 
change  only  the  cyclic  distribution  of  the  transient  reverse 
flow,  and  not  its  characteristic  development.   It  may  be 
generally  said,  that  of  the  three  prevalent  waveform  phe- 
nomena, only  transient  flow  reversal  significantly  changes 
the  character  of  the  instantaneous  velocity  profile. 
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C.   DEFINITION  OF  SEPARATION 

The  importance  of  separation  in  steady  flow  lies  in 
the  fact  that  its  occurrence,  as  characterized  by  the 
initial  disappearance  of  shear  at  the  wall,  marks  the 
beginning  of  a  process  that  ultimately  results  in  a 
streamline  displacing  region  of  wake  flow.   The  essential 
feature  of  this  process  is  the  development  of  upstream 
fluid  convection,  or  reverse  flow,  in  the  boundary  layer, 
resulting  in  the  production  of  eddies  which  induce  wake 
formation.   To  be  relevant,  any  process  or  phenomenon 
defined  as  separation  in  oscillating  flow,  must  be  sim- 
ilarly related  to  the  formation  of  a  region  of  wake  flow. 

The  findings  of  this  investigation  amply  verify 
previous  analytical  predictions  that  the  steady  flow 
definition  of  separation  is  inapplicable  in  oscillating 
flows.   Reverse  flow,  and  consequently  zero  wall  shear 
profiles,  were  observed  to  occur  periodically  at  all 
points  on  the  surface  in  the  adverse  pressure  gradient 
regime,  upstream  of  the  region  of  wake  formation.   These 
transient  flow  reversals  caused  no  observed  perturbations 
either  in  the  boundary  layer  or  the  freestream.   Further, 
it  can  reasonably  be  assumed  that  transient  flow  reversal 
occurs  even  in  the  zero  pressure  gradient  regime  for 
many  operating  conditions.   This  phenomenon  clearly  pre- 
cludes the  use  of  the  steady  flow  definition,  since  no 
initial  occurrence  of  a  zero  wall  shear  profile  or  flow 
reversal  can  be  logically  related  to  wake  formation. 
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Since  time  dependency  is  the  complicating  feature  in 
oscillating  flow,  the  use  of  a  time-average  or  mean  pro- 
file in  defining  separation  seems  logical.   However,  as 
has  previously  been  discussed,  mean  profiles  obtained  with 
a  hot-wire  anemometer  are  erroneous  when  transient  flow 
reversals  are  present. 

After  considering  all  such  possible  alternatives, 
only  one  observable  flow  phenomenon  is  clearly  related 
to  the  development  of  wake  flow.   This  is  the  initial 
occurrence  of  zero  velocity  or  reverse  flow  instantaneous 
profiles  throughout  the  entire  cycle  of  oscillation, 
or  "continuous  flow  reversal."   For  all  operating  con- 
ditions observed,  the  region  of  wake  formation  invariably 
began  in  the  immediate  neighborhood  of  the  initial  point 
of  continuous  flow  reversal.   Moreover,  it  was  found  that 
this  point  could  be  rigorously  and  uniquely  located  for 
all  observed  operating  conditions. 

The  following  definition  is  therefore  proposed: 

Laminar  boundary  layer  separation  in  oscillating  flows 
is  defined  as  commencing  with  the  initial  occurrence  of 
zero  velocity  or  reverse  flow  at  some  point  in  the  velocity 
profile  throughout  the  entire  cycle  of  oscillation. 

This  definition  has  the  following  important  features: 

1.   The  phenomenon  described  has  been  experimentally 
observed  to  be  directly  related  to  the  produc- 
tion of  wake  flow. 
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2.  Separation,  thus  defined,  occurs  at  a  unique 
point  in  the  flow  field. 

3.  The  definition  may  be  extended  to  steady  flow 
with  no  ambiguity. 

4.  The  definition  is  useful  in  analytic  as  well 
as  experimental  investigation. 

The  separation  point  then  is  defined  as  that  point  where 
the  first  indications  of  time-wise  continuous  boundary 
layer  flow  reversal  are  observed. 

A  suitable  criterion  for  the  recognition  of  separation 
is  the  observation  of  a  velocity  gradient  at  the  wall  which 
is  less  than  or  equal  to  zero  throughout  the  entire  cycle 
of  oscillation, 


~£  ■  o   for  y  =  0   and   0  <  t/m   <    1 
dy  —  1 


Figure  28  illustrates  typical  cyclic  variation  of  the 
instantaneous  velocity  profile  at  separation.   Figure  29 
presents  the  associated  waveform  photographs  from  which 
the  profiles  were  derived.   Figure  30  shows  the  change  in 
character  of  the  instantaneous  velocity  profile  sequence 
at  separation  accompanied  by  the  presence  of  strong  flow 
reversal.   Figure  31  presents  the  associated  waveforms. 

Figures  32  through  59  show  instantaneous  profile 
envelopes  at  separation  for  a  wide  range  of  operating 
conditions.   Figure  50  illustrates  the  influence  of  large 
phase  leads  on  the  profile  envelope.   These  figures  are 
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representative  of  27  sets  of  separation  point  data  obtained 
from  the  original  82  profile  sequences.   Once  the  separation 
criteria  were  established  from  this  information,  the 
pressure  distribution  on  the  model  was  varied  and  an  addi- 
tional 36  separation  points  located  by  direct  reading  of 
the  data  displays.   Steady  flow  separation  points  were 
located  in  the  usual  manner  for  each  of  the  combinations 
of  mean  velocity  and  pressure  distribution  employed  in 
the  investigations.   These  are  listed  in  Table  II.   The 
change  in  separation  point  location  due  to  the  imposition 
of  oscillation  was  then  determined.   All  oscillating  flow 
separation  data  are  listed  in  Table  III. 

D.   CORRELATION  OF  SEPARATION  POINT  ALTERATION 

The  correlation  between  oscillating  flow  separation 

point  location  and  flow  parameter  variation  has  been 

obtained  by  a  "least  squares"  fit,  the  details  of  which 

are  outlined  in  Appendix  B. 

In  terms  of  dimensionless  deviation  from  steady  flow 

separation,  the  following  dimensionless  correlation  was 

obtained: 

a      ss    s    ,     ,-.  ,  A7  .T  -1.77   -.24  -.26  --.035 

A   =  =  1.13  x  10   N_       N_  N-     N,, 

s      x  Re       P  f      A 
ss 

The  correlation  represents  the  data  to  within  three  per 
cent  and  the  average  scatter  is  less  than  one  per  cent. 
This  indicates  a  regularity  of  occurrence  that  strongly  re- 
inforces the  validity  of  the  present  definition  of  separation. 
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TABLE  II 
STEADY  FLOW  SEPARATION  DATA 


U       (ft/sec)  N„    xl0~5  Nn  x      /c 

oo  /  Re  p  ss/ 


Pressure  Distribution  I 


15 
20 
25 


2.427 

.3016 

.724 

3.033 

.2891 

.679 

3.632 

.1083 

.650 

Pressure  Distribution  II 


15 
20 
25 


2.163 

.4485 

.645 

2.746 

.3362 

.614 

3.339 

.2501 

.598 

Pressure  Distribution  III 


15 
20 
25 


2.602 

.3573 

.776 

3.267 

.2478 

.731 

3.885 

.0497 

.695 
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The  cumulative  effects  of  possible  experimental 
inaccuracies  on  the  correlation  function  were  considered. 
The  per  cent  uncertainty  inherent  in  the  correlation 
function  was  estimated  using  the  method  of  Kline  and 
McClintock  [34].   The  following  is  a  brief  outline  of 
this  procedure: 

^00 

Setting  K  =  N   — ,  and  rearranging  the  correlation 

o 
function; 

r=n -1-77  V24  n--26  n;035 

Re       P      f      A 

-1.25   -1.77   -.26   1.51   -.24    -.24    .24  ^.035  _-.035 
=  (UJ      (x)      (u>)     (v)     (K)    (xq)    (qj    (U)    (U) 

Differentiating  logarithmically  and  dividing  by  R; 

M=-1.25  ^-1.77  M  -  .26  ^+  1.51  *±  - -.035  B 

R  n  x         oj  V  U 

oo 

Taking  the  root  square  (the  most  probable  combinatorial 
sum)  ; 

as  =  [(-1.25?^)2+(-i.77aV+<-.26aV+ +(-.0  35aj)2]1's 

R  U  x  °°  U 

00 

The  respective  uncertainty  intervals  were  estimated  to 
be; 

00  7\\)  ^oo 

=  .02  —  =  .01  -=:—  =  .05 


3X  =  .002  £-  .05  22|-  .05 


x     *~~~  K 


U 
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^=    .001 

0) 

9x 

o 

X 

o 

Phus  ; 

3R  _ 
R 

.034 

.002  ^-   =  .02 

U 


Hence,  there  is  an  estimated  uncertainty  of  3.4  per 
cent  in  the  correlation  function's  representation  of  the 
actual  physical  phenomena  within  the  parametric  limits  of 
this  investigation,  namely; 

1.950xl05  <  N'   <  3.637xl05 
—  Re  — 

1.852xl0"5  <  Nf  <  15.709xl0-5 
.0100  <  Np  <.  .2723 
.03  <  NA  <  .44 

For  this  range  of  flow  parameters  the  correlation 
function  indicates  that  the  oscillating  flow  separation 
point  moves  toward  the  steady  flow  value  with  increasing 
frequency  and  decreasing  amplitude,  as  one  would  intui- 
tively expect.   It  is  apparent,  however,  that  within  this 
parametric  range,  the  effect  of  amplitude  variation  on 
separation  point  location  is  less  than  the  estimated  degree 
of  uncertainty. 

Figures  59  through  62  are  parametric  plots  of  the 
correlation  function.   Dashed  lines  indicate  extrapolations 
beyond  the  range  of  observed  data. 
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FIGURE  59 

SERRATION   POINT  ALTERATION 
AS  A  FUNCTION    OF    PRESSURE     PARAMETER 
AND    REYNOLDS    NUMBER  WITH 
Nf     =  I.Oxl65 
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FIGURE  60 

SEPARATION    POINT    ALTERATION 
AS  A   FUNCTION    OF    PRESSURE     PARAMETER 
AND    REYNOLDS   NUMBER  WITH 
Nf  =  5.0  x  I05 
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FIGURE  61 

SEPARATION    POINT  ALTERATION 
AS  A    FUNCTION  OF     PRESSURE  PARAMETER 
AND    REYNOLDS  NUMBER  WITH 
Nf  *IO.Oxl<r 
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FIGURE  62 

SEPARATION    POINT  ALTERATION 
AS  A   FUNCTION   OF  PRESSURE    PARAMETER 

AND    REYNOLDS  NUM^R  WITH 
N    =I5.0xI0"5 
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VI.   DISCUSSION 

A.   PERIODIC  WAVEFORM  DISTORTION 

During  the  course  of  investigation,  several  general 
trends  in  the  occurrence  of  the  two  previously  described 
periodic  waveform  effects  were  observed. 

1 .   Amplification 

Maximum  amplification  at  any  chordwise  position, 
for  fixed  operating  conditions,  generally  occurred  between 
the  20  per  cent  and  50  per  cent  freestream  velocity 
laminae  in  the  mean  profile,  similar  to  the  findings  of 
Hill  [18]  and  Hori  [19].   With  all  other  flow  parameters 
fixed,  maximum  amplification  varied  inversely  with  the 
amplitude  of  the  freestream  oscillation.   Maximum  amplifi- 
cations of  40  per  cent,  100  per  cent,  and  200  per  cent 
were  attained  for  values  of  NA  in  the  range  of  .30,  .10, 
and  .03  respectively.   In  all  instances  of  its  occurrence, 
amplification  was  observed  to  increase  in  magnitude  with 
increasing  streamwise  coordinate,  but  exhibited  no  identi- 
fiable trend  with  mean  velocity  or  pressure  gradient. 
Figure  6  3  presents  values  of  maximum  amplification  as  a 
function  of  freestream  amplitude.   Comparison  with  Hill's 
findings  indicates  that  in  the  present  work,  either  the 
larger  adverse  pressure  gradient  or  higher  frequencies,  or 
both,  tend  to  induce  greater  amplification  in  the  case  of 
small  amplitudes.   At  larger  amplitudes,  the  present 
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values  appear  to  be  in  closer  agreement  with  Hill's  zero 
pressure  gradient  results. 
2.   Phase  Shifts 

The  phase  shifts  observed  were  invariably  phase 
leads,  as  originally  predicted  by  Lighthill  [4],  and 
observed  by  Nickerson  [7]  and  Hill  [18]  in  oscillating 
Blasius  flow.   For  values  of  N   of  the  order  of  .30,  the 
phase  lead  was  almost  non-existent,  and  reached  a  maximum 
of  only  10  degrees  when  it  did  occur.   Phase  leads  of  up 
to  30  degrees  were  observed  for  values  of  N   of  approxi- 
mately  .10.   For  values  of  N   of  the  order  of  .03,  the 
maximum  observed  phase  lead  was  approximately  45  degrees 
for  most  operating  conditions.   A  significant  exception  was 
noted  at  these  smallest  amplitudes;  at  a  Reynolds  Number 
of  .98xl05  per  foot  and  Nf  of  5.144xl0~5,  phase  leads  of 
up  to  180  degrees  were  repeatedly  observed,  accompanied  by 
amplifications  of  up  to  200  per  cent.   Nothing  approaching 
effects  of  this  magnitude  was  observed  for  any  other 
operating  condition.   Similar  to  amplification,  phase  leads 
exhibited  no  dependency  on  mean  velocity  or  frequency, 
other  than  in  the  exceptional  case  mentioned  above.   Figure 
64  presents  values  of  maximum  phase  lead  as  a  function  of 
freestream  amplitude.   Comparison  with  zero  pressure 
gradient  results  indicate  that  with  the  exception  of  the 
case  mentioned  above,  the  presence  of  the  adverse  pressure 
gradient  tends  to  suppress  phase  leads  for  freestream 
amplitudes  greater  than  .10. 
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B.   TRANSIENT  FLOW  REVERSAL 

It  has  been  generally  observed,  in  the  case  of  steady 
flow,  that  boundary  layer  flow  reversal  causes  vortex 
formation,  which  in  turn  induces  wake  flow  [32].   Observing 
the  phenomena  of  transient  flow  reversal,  one  quite 
naturally  questions  the  absence  of  disturbances  in  the 
boundary  layer  due  to  the  reversal.   One  physically 
reasonable  explanation  is  that  the  transient  reversals 
are  of  such  short  temporal  duration,  that  local  vortices 
do  not  have  time  enough  to  grow  to  sufficient  strength  to 
disturb  the  flow.   It  has  been  observed  that  in  the  case 
of  boundary  layer  separation  on  accelerating  bodies,  a 
finite  time  increment  exists  between  initial  flow  reversal 
and  vortex  formation  [32].   This  line  of  reasoning  also 
suggests  that  at  very  low  frequencies,  with  relatively 
long  periods  of  "transient"  flow  reversal,  wake  formation 
might  conceivably  begin  with  the  initial  occurrence  of 
reverse  flow.   Such  a  possibility  is,  of  course,  compatible 
with  the  definition  and  theory  of  "quasi-steady"  oscillat- 
ing boundary  layer  flows  advanced  by  Lighthill  [4],  Lin  [5] 
and  others.   This  picture  is,  however,  not  in  consonance 
with  Stuart's  [28]  predictions  of  transient  flow  reversal 
in  a  zero  pressure  gradient,  which  the  results  of  this 
investigation  strongly  suggest  actually  occur. 

Following  the  previous  line  of  reasoning,  this  implies 
the  possibility  of  separation  occurring  on  a  flat  plate 
in  low  frequency,  large  amplitude  oscillating  flow,  a 
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violation  of  the  equations  of  motion.   It  seems  obvious 
that  much  more  knowledge  of  the  interaction  of  time- 
dependent  momentum,  pressure,  and  viscous  forces  is  required 
before  phenomena  such  as  transient  flow  reversal  can  be 
fully  understood 
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VII.   CONCLUSIONS 

From  the  results  obtained,  the  following  conclusions 
may  be  drawn: 

1.  The  first  indications  of  wake  formation  occur  at 
that  streamwise  point  in  the  flow  where  the  normal 
velocity  gradient  at  the  wall  is  observed,  for  the 
first  time,  to  be  less  than  or  equal  to  zero 
throughout  the  entire  cycle  of  oscillation. 

This  point  may  then  be  logically  defined  as  the 
"separation  point,"  marking  the  beginning  of 
boundary  layer  separation. 

2.  The  presence  of  an  oscillating  velocity  component, 
superimposed  on  the  mean  flow,  causes  separation 
to  occur  upstream  of  the  steady  flow  separation 
point. 

3.  Increasing  frequency  of  oscillation  tends  to  cause 
the  separation  point  to  move  downstream  toward  the 
steady  flow  location. 

4.  For  the  range  of  flow  parameters  observed,  the 
effects  of  amplitude  variation  on  separation  point 
location  are  less  than  the  estimated  experimental 
uncertainty. 

5.  Transient  boundary  layer  flow  reversal  occurs 
throughout  the  adverse  pressure  gradient  regime 
for  amplitude  parameters  of  the  order  of  .10  or 
greater,  without  inducing  wake  formation. 
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APPENDIX  B 
MATHEMATICAL  OUTLINE  OF  PARAMETRIC  CORRELATION 

The  correlation  between  oscillating  flow  separation 
point  location  and  flow  parameter  variation  was  obtained  in 
the  following  manner: 

It  was  desired  to  determine  an  expression  of  the  form: 

A   =  _ss S  *  K  N  a   b    c   d  (1) 

s       x  Re   P    f    A  ' 

ss 

4s  =  (As  +  E)  =  K  NR*  Npb  NfC  N/  (la) 

where  A  ,  N   ,  N  ,  N->  and  N   are  experimentally 

obtained  and  K,  a,  b,  c,  and  d  are  constants  determined 

in  such  a  way  as  to  minimize  the  error  (e)  in  the 

correlation. 

Simplifying  the  notation: 


NA  =  ^5 
a   =  a2 

b   =  a. 


A 
s 

" 

X 

^, 

~ 

A 
s 

s 

X 

K 

= 

al 
e      = 

al 
*1 

NRe 

= 

y2 

NP 

= 

^3 

Nf 

= 

^4 

c   =  a4 


d   =  a5 
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equation  (la)  becomes: 

al   a2   a3   a4   a5 
x  =  y±     y2  y3  y4  y5 

or,  in  logarithmic  form: 


£n  x  =  Z  a.  in  y .  i  =  1,2,... ,5     (2) 

i   -1      1 

now  considering  all  the  separation  point  data: 


Z  Jin  x  .  =  Z  ( Z  a .  £n  y .  . )      i=l,2,...,5     (3 

D     J         D   i  J 

j  =  1,2,.. .,63 


This  expression  was  noted  to  be  adaptable  to  treatment 
by  application  of  a  modified  "least  squares"  curve  fit 
method.   Accordingly,  "residuals,"  in  terms  of 
observed  separation  point  variation  were  defined: 


Z  R.  =  i-  Z  (Z  a.  in   y.  .  -  Jin  x.)2  (4) 

j  ]   2  j  i  i   ^i]      : 


minimizing  the  residual: 

~-  (Z  R.)  =  Z(Z  a.  iny..-in   x .  )  i      y,   =  0 
9ak  \   :    j  i   X      D       D   n    : 

xC  =   J-,Z,...,!Z) 


or: 


Z  Z  a.  in   y.   £n  v.   -  Z  in   x .  in  y.   ,  =  0       (5) 
31         J  j 
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in  matrix  notation,  with: 

{A}  =  a  5x1  matrix  of  "a." 

[Y]  =  a  5x63  matrix  of  "Jin  y  .  " 

{X}  =  a  6  3x1  matrix  of  "In   x." 

: 

equation  (5)  becomes: 

[Y] [Y]T{A}  -  [Y]{X}  =  0  (6) 

(A)  =  [YY1]   {YX}  (7) 


equation  (2)  becomes 


{X}  =  [Y]T{A}  (8) 


and  from  (la) : 

{X}  -  {X}  =  {e}  (9) 

Using  standard  matrix  algebra  subroutines,  a  computer 
solution  of  equation  (7)  was  obtained,  resulting  in: 

al  7 

K  =  e  -1  =  1. 13x10 ' 

a  =  -1.77 

b  =  -.24 

c  =  - . 26 

d  =  .035 

The  desired  correlation  function  was  then: 


a  =  !lsl^=i.i3xio7  „  -i-"  n--24  n:-26  n:035 

s     x  Re       P      f      A 

ss 
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Using  the  solution  of  equation  (7) ,  equations  (8) 
and  (9)  were  successively  solved,  resulting  in  the 
determination  that  the  above  function  represents  the 
data  to  within  three  per  cent,  with  average  scatter 
of  only  one  per  cent. 
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